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QO ■ Abstract 

' The effects of intermediate charmed mesons on charmonium transitions with the emis- 

sion of one pion or eta are studied systematically. Based on a non-relativistic effective 
field theory we show that charmed meson loops are enhanced compared to the correspond- 
ing tree-level contributions for transitions between two 5*-wave charmonia as well as for 
' transitions between two P-wave charmonia. On the contrary, for the transitions between 

one S'-wave and one P-wave charmonium state, the loops need to be analyzed case by 
case and often appear to be suppressed. The relation to and possible implications for an 
effective Lagrangian approach are also discussed. This study at the same time provides a 
cross check for the numerical evaluations. 
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1 Introduction 



Since the discovery of the J/tp more than thirty years ago, the decays of heavy quarkonia 
have played an important role in the physics of quarks and hadrons. During the past decades, 
experimental studies of the heavy quarkonia at CLEO, Belle, BaBar, CDF, DO, and BES- 
II have provided great opportunities for examining many interesting properties of Quantum 
Chromodynamics (QCD). At the present stage, BES-III jJJ has accumulated the largest data 
samples for J/t/j and ip' decays, and PANDA [2 J plans to accumulate data for charmonia which 
cannot be produced directly in electron-positron annihilations. These facilities will deepen our 
understanding of the charmonium physics, and hence the non-perturbative aspects of QCD. 
Although many theoretical investigations have been performed in the past thirty years (for 
comprehensive reviews, see Refs. [3j [JJ 2]), there remain many mysteries in charmonium 
physics to be settled. On the contrary, due to the new experimental data with unprecedented 
statistics, many new interesting problems have appeared, e.g. the nature of many of the 
new X, Y, Z resonances discovered in the charmonium mass region has still not been well 
understood (for recent reviews, see, e.g. Refs. El [7J El [9] ) . 

Furthermore, various recent phenomenological calculations suggest that charmed meson 
loops may play an important role in the decays of heavy quarkonia (for an overview, see 
|10|). For instance, using an effective Lagrangian approach (ELA), intermediate heavy-meson 
loop contributions are found to be essential for understanding the puzzling ^(3770) non-DD 
decays [HI [12] . They are also important in the J / ip decays into a vector and a pseudoscalar 
meson |13j and in the Ml radiative transitions between two charmonia [13] • Besides, using the 
on-shell approximation, the bottom meson loops were suggested to make the T(5S I ) transitions 
to the lower T states with the emission of two pions |15] or one r] [16] different from those 
of the T(4S). The inclusion of intermediate heavy mesons in heavy quarkonium transitions, 
sometimes called coupled-channel effects, has been noticed for more than twenty years [17} 
ITS], 119]. Also, the effect of the mass differences between the neutral and charged mesons in the 
intermediate states (i.e. in the meson loops) plays a role in other isospin breaking processes. 
This effect, known to be of particular importance near the continuum thresholds, was already 
studied in the rf decays [20], and in the decays (p -> wtt° [HI [22], J/ip -> (prjir (22 [23], and 
L> s0 (2317) -> D s ir° [23 EH [27] . 

An often used formalism dealing with the hadronic transitions between two heavy quarko- 
nia is the QCD multipole expansion (QCDME) [2SJ [291 EQ] . The QCDME is based on the 
assumption that the emitted gluons are soft so that their wavelengths are much larger than 
the size of a heavy quarkonium. As a result, a multipole expansion similar to that in classical 
electrodynamics can be performed. The soft gluons then hadronize into light meson(s), for in- 
stance the pion(s) or eta, and the matrix elements may be worked out using soft pion theorems. 
A schematic diagram for the multipole transition from a heavy quarkonium to another one 
with the emission of one pion is plotted in Fig. [TJa). However, this ansatz clearly misses the 
contribution from intermediate mesons. This can be understood as a heavy quarkonium can 
couple to a heavy meson and heavy anti-meson pair through the non-perturbative production 
of a light quark and anti-quark pair, see Fig. [JJb). In Refs. |19| 13 lj . within the framework of 
QCDME, the intermediate heavy meson effects were considered to account for non-multipole 
effects in the di-pion transitions between two charmonia, the ip states, or bottomonia, the T 
states. Better agreement with the experimental data was obtained. 

In light of these phenomenological indications, it is important to have a theoretical for- 
malism which has a controlled uncertainty to study the effects of heavy meson loops in the 
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Figure 1: Schematic diagrams of the QCDME mechanism (a) and the non-multipole (b) effects 
of the intermediate heavy meson loops for heavy quarkonium transition with the emission of 
one pion. 



transitions of heavy quarkonia. In Ref. [32J, such a non-relativistic effective field theory 
(NREFT) formalism was constructed to investigate the charmed meson loops in the decays 
tp' — > J/ip7T°(r]). Because the i/j' and J ftp are isospin and SU(3) flavor singlets, the decay 
process ip' — > J/ipir° violates isospin symmetry and tp' — > J/tpr] violates SU(3) symmetry. 
Isospin symmetry can be violated by both electromagnetic (e.m.) effects and the mass dif- 
ference between the u and d quarks. Because the e.m. effects are small |33|, [M] (it can 
be easily shown in the framework of chiral perturbation theory (CHPT) with virtual pho- 
tons [351 EE] , see Section 12. 2p , these two decays were used to extract the light quark mass 
ratio m u /m<i [371 ESI EHl HOI HI]. The relation between the ratio of the decay widths 
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where F n ^ and M w r^ are the decay constant and mass of the pion (eta), respectively. The 
extracted quark mass ratio using Eq. (pQ) and the recent measurements of the decay widths 
from the CLEO Collaboration [32], the BES Collaboration [13] and the Particle Data Group 
(PDG) fit |44| are listed in Tabled! Comparing with the result obtained using the Goldstone 
boson masses from leading order (LO) CHPT 



M 2 K+ 



M 2 R0 + 2M 2 



M 1 



M K0 



M 2 K+ +M^ 



0.56, 



(2) 



the discrepancy is striking. We remark that there might be sizable higher order corrections 
to this LO result. The up-to-date knowledge of the light quark mass ratio from various 
determinations including lattice calculations (but excluding ijj' decays) was summarized by 
Leutwyler as m u /md = 0.47 ± 0.08 |47] , The relatively large uncertainty given here thus 
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RjP/h rn u /m d 


CLEO 02] 
BES [43] 
PDG fit [S] 


(3.88 ±0.23 ±0.05)% 0.40 ± 0.01 
(4.8 ±0.5)% 0.35 ±0.02 
(4.0 ±0.3)% 0.39 ±0.02 



Table 1: The light quark mass ratio m u /md extracted using Eq. ([I]) from the recent experi- 
mental measurements by different collaborations. 



provides an overlap with the results quoted in Table HJ however, only at the very low end. 
In Ref. |32] it was stressed that for the mentioned transitions the effects of charmed meson 
loops, ignored in the previous analyses, should be sizable. The charmonia ^ and J/ip couple 
to charmed and anti-charmed mesons, and the pion is emitted from one intermediate charmed 
meson. The proper expansion parameter is the velocity of the intermediate meson, which for 
below threshold decays is defined via the analytic continuation of the standard definition, 
namely v = yj —E /Mp, with E measured relative to the open charm threshold. We find 
v 0.5 for most of the decays studied. It is found that loops are enhanced by a factor 
of 1/v compared to the tree-level contribution where the pion is emitted directly from the 
charmonium. Therefore the dominant (LO) contributions to these decays come from the loops 
instead of from the tree graphs which are proportional to the quark mass differences directly, 
and hence the extraction of quark mass differences from ip 1 decays mentioned above is not 
reliable. Stated differently: quarkonium decays could only be used to extract the light quark 
mass ratio, if either the loop contributions could be controlled quantitatively (at present their 
uncertainty is quite sizable — see discussion below), or if loop contributions are suppressed. 
It turned out that the enhancement of the loops in case of the tp' to J/ tjj transitions emerges 
only because the transitions at hand violate isospin or SU(3) symmetry. The power counting 
is discussed in detail in Sec. 13.41 

In Ref. [48], the same NREFT is applied to the decays tp' — > h c 7T° and r( c — > XcO^ - 
As a consequence of the quantum numbers of the charmonia involved, in these two decays 
the loop contributions are highly suppressed, and hence the tree-level terms, i.e. the quark 
mass difference terms, dominate the decay amplitudes. Unfortunately there is no charmonium 
transition, where one can exploit this observation in order to extract the light quark mass ratio, 
since typically the phase space available is insufficient for an r\ in the final state. However, 
in the bottomonium system analogous transitions appear to exist [49J, and will allow for 
the mentioned analysis. This illustrates that the effective field theory at hand predicts a 
highly non-trivial pattern for the loop contributions in different decays that can be tested 
exp er iment ally. 

In this paper, we will systematically investigate the charmed meson loop contributions 
to the transitions between two charmonia with the emission of one light pseudoscalar me- 
son. We restrict the charmonia to S and P-wave states with radial quantum number n less 
than or equal to 2. In other words, we will consider the transitions between or within the 
following charmonia spin multiplets: { J/V', 77c}, {Xe0> Xd, Xc2, h c } with n = 1, and {ip',rj' c }, 
{Xco> Xci> Xc2' ^4} with n = 2. Charge conjugation allows for the emission of one light pseu- 
doscalar meson between two charmonia with the same value of C. Considering further the 
constraints from parity conservation, all the allowed transitions are plotted in Fig. [21 and 

In the figure, the masses of the so far unobserved (or unidentified) 2P charmonia are taken from a quark 
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the following will be considered in the paper: 

1) Transitions between two S- wave charmonia: ip' — > J/ipir°, and ip' — > J/iprj. 

2) Transitions between one S- wave and one P-wave charmonium: ip' — > h c 7r°, h c — > J/ipn , 
K ->• ip'n , Vc -> Xco^r , Xco -> VcK°, and x' c o -> r ]'c 7T °- 

3) Transitions between two P-wave charmonia: x' c o ~^ Xci^°, x' c i ~^ Xcj^° {J = 0,1,2), 
Xc2 -> Xci(2)K°, and ^ -> /i c vr°. 

In fact, as plotted in Fig. [2J the decays /i' c — > J/ipir°(r]) and x' c o — ^ W 71 " ^) can also occur. 
However, the mass difference between the initial and final charmonium exceeds 800 MeV, of 
order 0(A X ), with A x ~ 1 GeV denoting the typical hadronic scale. Since the chiral expansion 
is an expansion in p/A x , with p denoting a typical momentum or mass, for those energies 
the chiral expansion is not expected to converge any more. We therefore do not consider 
these transitions. There could also be .D-wave transitions Xc2 r ?c 7 i"°( ? ?) an d x' C 2 ~^ r ]'c 7T ° ( n °t 
shown in the figure). However, their partial decay widths would be too small to be detected 
in the near future because of the D-wave suppression, the isospin or SU(3) breaking, and 
small phase space. They will also not be considered here. 

The paper is organized as follows. In Section [21 the tree-level chiral effective Lagrangians 
and the resulting amplitudes are given for all the decays discussed in the paper. The isospin 
and SU(3) breaking are given by quark mass differences. Consistent with earlier analyses, 
e.m. contributions are found to be small and can be neglected. Various aspects of the 
charmed meson loops will be discussed in Section [3J Section H] is devoted to the results of 
the meson loops for various transitions in the NREFT. In Section [5j a detailed parallel study 
of the meson loop transitions in the framework of the ELA is presented and the results are 
compared to those from the NREFT. A brief summary is given in the last section. Various 
technicalities such as the utilized loop functions, the decay amplitudes from the meson loops 
and the ingredients of the ELA are given in the Appendices. 

2 Effective Lagrangians for tree-level diagrams 

Isospin breaking has two sources. One is the mass difference between the up and down quarks, 
and the other one is of e.m. origin because photons do not have definite isospin. This section 
is devoted to the construction of the LO chiral Lagrangians for the tree-level diagrams of the 
transitions considered in this paper. Both the quark mass difference and the e.m. effects will 
be taken into account. 

The chiral effective Lagrangians are of the most general form which is invariant under 
the transformations of SU(3)ixSU(3)^, parity and charge conjugation. The charmonia are 
treated as matter fields, and the pion and eta are the Goldstone bosons of the spontaneous 
breaking of SU(3)lxSU(3)# down to its vector subgroup SU(3)y. The charmonia are SU(3) 
singlets, so they do not change under the chiral transformation. The Goldstone boson fields 



model calculation considering the color-screening effect due to the light quark and anti-quark pair creation [50] . 
Note that these values are only used for illustration. 
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Figure 2: All possible 5 and P-wave transitions among the ground state and the first radial 
excited S and P-wave charmonia with the emission of one light pseudoscalar meson (77° , 77) . 
The unobserved resonances and decay modes are plotted in dashed. The masses for the 
unobserved resonances are taken from predictions in Ref. [50 1 . The thresholds for the DD, 
DD* and D*D* are represented by the dotted horizontal lines. The solid lines represent the 
measured decays. The dashed and dotted lines are not yet measured with only the former 
discussed in the paper and the latter are beyond the range of the applicability of the chiral 
EFT (but can be considered in the model-dependent ELA also discussed here). Thick lines 
indicate transitions with enhanced charmed meson loops. 



where we have approximated the 77 as one element of the octet SU(3) representation, are 
collected in u = exp (i(j)/V2F) with F being the pion decay constant in the chiral limit. 
Under the transformation of SU(3)lxSU(3)^, we have 



u 



Ruh ] = huh\ (4) 



where h is the compensator field. It is convenient to construct chiral Lagrangians using 
operators whose chiral transformation is O — > hOh^ . The following such building blocks will 
be used 



X± = u] XU ] ± ux'u, 

Q± = \ ( ui Q u ± u Q u] ) . ( 5 ) 

where the diagonal quark mass matrix and the charge matrix are 

X = 2B ■ diag (m u , m d , m s ) , 

Q = e • diag (2/3, -1/3, -1/3) , (6) 
in terms of Bq = \(0\qq\0)\/F 2 and the elementary electric charge e (e > 0). 
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In the heavy quark limit with txiq —> oo, the coupling of a heavy quark to a gluon is 
spin-independent [51]. As a result, there is a spin symmetry in that limit, and the heavy 
quarkonia, which differ from each other only in the total spin of the heavy quark and anti- 
quark, can be grouped into the same spin multiplet. It is then convenient to introduce a single 
field for a spin multiplet of heavy quarkonia |52l Iq"5] using the trace formalism proposed for 
single- heavy mesons [54, 55J. In this way, the consequence of the heavy quark spin symmetry 
can be obtained automatically by evaluating a trace in spinor space. The construction of 
charmonium fields with arbitrary orbital angular momentum I in the trace formalism can 
be found in Ref. [53J. Since we are dealing with the transitions between two charmonia, 
the heavy quark four-velocity is conserved up to higher order corrections. In this case, it is 
convenient to use the two-component notation as introduced in Ref. [56J. For doing that, the 
four- velocity is chosen to be = (1,0). In the two-component notation, the field for the 
5-wave charmonia reads 

J = ip ■ a + rjc, (7) 
with if) and rj c annihilating the J/tp and r] c states. The field for the P-wave charmonia isH 

X* = <r J i-X% ~ ^ ljk Xci + 771^X00) + K, (8) 

where x&i Xci> XcO and h c annihilate the Xc2, Xci, XcO and h c states, respectively. x % di is a 
symmetric and traceless tensor. 

The quantum numbers of the charmonia determine their parity and charge conjugation 
transformation properties. The parity transformations for the charmonia fields are given by 

J 4 -J, x'^x\ (9) 
and the charge conjugation transformations are given by 

C T ~* 

J — > a 2 J 01 = — ip ■ cr + T] c , 

X 1 4 -<7 2X V 2 = j [~X% ~ ^ jk Xci + ^Xco) - K, (10) 

where J T is the transpose of J. Denoting the rotation in the SU(2) spin space of the heavy 
quark (anti-quark) by 5 (S), the transformation of the charmonium fields reads 

J 4 SJS^ x i ^ s^sl (11) 

The transformation for a heavy quarkonium field with arbitrary orbital angular momentum 
is given in [53] in four-component notation. In two-component notation, the transformation 
properties for the Xcj fields can be found in Ref. [57] . 

2.1 Quark mass difference 

In the transitions considered in this paper, one pion or eta is emitted. Therefore we need to 
construct the chiral Lagrangian using the external field X- > which is proportional to the light 

2 The sign convention for the Xc2 and Xrf fields are different from those in Ref. [57] . 



7 



quark mass matrix and contains an odd number of the Goldstone bosons. Under parity and 
charge conjugation, X- transforms as 



V 




(12) 



X- -> -X- = XL, 



respectively. 

The LO Lagrangian for the transitions between two S-wave charmonia can be constructed 
considering parity conservation, which requires the presence of a derivative, charge conjuga- 
tion, chiral symmetry and Galilean invariance, 



where (...) denotes the trace in spinor space, the subscript a = u, d, s is a flavor index, and aa 
as a sum over it denotes the trace in flavor space. Similarly, the Lagrangian for the transitions 
between one S- and one P-wave charmonium states is 



These Lagrangians were first proposed in Ref. [58] in four-component notation. Note that due 
to the presence of a Pauli matrix between the two heavy quarkonium fields in Eqs. (|13ll4p . 
the heavy quark spin symmetry is violated. On the contrary, the Lagrangian Cpp preserves 
the spin symmetry. 

The goal of the present work is to set up an effective field theory that allows one to 
systematically study both loop as well as tree level transition amplitudes. To prepare for 
this we need to assign an order of magnitude estimate to the coupling constants given above. 
They may be determined in principle as the result of some matching procedure between the 
hadronic matrix elements and the more fundamental quark-gluon dynamics calculated within 
(potential) nonrelativistic QCD ((p)NRQCD) — for recent reviews, see Ref. [59J, El [H H] - For 
instance, the coupling constant A in the Lagrangian Eq. (|13p has a mass dimension —2. 
There are several different scales in the physics related to heavy quarkonia. They are the 
heavy quark mass rriQ, the momentum uiqvq, the inverse of which sets the length scale of a 
heavy quarkonium, and the energy scale vhqVq [60], where vq denotes the velocity of the heavy 
quark within a heavy quarkonium to be distinguished from the velocity v of the heavy mesons 
in the loops to be introduced in the next section — for an estimate of the values of various 
scales in heavy quarkonia, one may refer to [ST]. In addition, there is the nonperturbative 
QCD scale Aqcd- I n this paper, we are considering the low- lying heavy quarkonia. For these 
states, it is believed that mgyq > Aqcdj which defines the weak-coupling regime [59] — e.g. 
with v\ ~ 0.3 and m c = 1.5 GeV we find mQVq ~ 450 MeV. Since most of the coupling 
constants introduced in the Lagrangians are dimensionful, they should have certain scaling 
properties expressed by the above mentioned scales. The tree-level Lagrangian describes a 
process with the emission of soft gluons, which then hadronize into a pion or an eta. As 
mentioned in the Introduction, the applicable regime of our effective field theory is limited 
to the transitions with the pion (eta) energy much smaller than A x , i.e. E^r^ < 600 MeV. 




(13) 




(14) 



and that for the transitions between two P-wave charmonia is 




(15) 
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Hence the energy of the emitted gluons should also be < 600 MeV. Therefore, the proper 
dimensionful parameter that sets the scale for this nonperturbative process should be either 
m c v 2 , which is sometimes called ultrasoft, or Aqcd- As mentioned above, the charmonia 
considered in this paper are weakly-coupled, i.e. m c v 2 c > Aqcd- So conservatively, one may 
take Aqcd to set the soft scale. Furthermore, since the transition violates spin symmetry we 
have to put in a factor Aqcd / m c to finally get 

1 AqCD - 1 (16) 



A QC D2m c m c 2ml ' 

where a factor of l/(2m c ) was introduced to make A have the correct dimension. In addition, 
assigning the pion decay constant as F ~ Aqcd, and the quark condensate as |(0|<jg|0)| ~ 
A qcd> one has 

Bdu = -y( m d ~ m u ) ~ 6, (17) 

with 5 denoting the quark mass difference. So using the expressions given in Table [2j the 
tree-level amplitude for a pionic transition scales as 

Mf r f e ~ -Lg<5, (18) 

IIL C 

where l/(2m c ) has been canceled by the factor ^MiMf due to nonrelativistic normalization. 

The dimension of the coupling constant 7 in the Lagrangian for the transitions between 
two P-wave quarkonia is the same as that of A. However, for these transitions, the spin 
symmetry is preserved as can be seen from Eq. (|15p . So in the scaling of 7 the suppression 
factor Aqcd/ttIc should not be present. This is the only difference from that of A. The 
dimension of C in Eq. (|14p is higher than that of A or 7 by one unit. Therefore, analogous to 
Eq. (|18j) . the scaling of the tree-level amplitude for a pionic transition between two P-wave 
charmonia and that for a transition between one S- and one P-wave states should be given 
by 

A"C P e ~ i^— q5 > ~ 6 - (19) 

Aqcd 



2.2 Virtual photons 

The e.m. effects come from virtual photons exchanged in the processes. The inclusion of 
virtual photons has been first considered systematically for three-flavor CHPT in Ref. |36j . 

Since the photons are virtual, we need to consider operators with at least two powers 
of electric charge. At 0(m^a), with a = e 2 /(47r) ~ 1/137 the fine structure constant, one 
virtual photon is exchanged. There are three types of operators, and we will discuss them 
one by one: 

1) The virtual photon is exchanged between light quarks, which is the standard case in CHPT 
with virtual photons. One needs quadratic combinations of the spurions and Q_ which 
act on the light quarks. The parity and charge conjugation properties of the Q± are given 
by [62] 

Q± 4 ±Ql, Q± 4 ±Ql. (20) 
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There is only one light pseudoscalar meson in the final states of all the transitions con- 
sidered in the paper. Q + and Q_ contain even and odd number of the Goldstone fields, 
respectively, and their expansion reads 



Q + = Q + O (0 2 ) , 

Q_ = *(Q0-^Q)+O(03). 



(21) 



Therefore, at 0(a), the possible virtual photon operators for the one pion (eta) emission 
transitions between charmonia are (Q+Q-) a a and (Q+) a a(Q-)aa- The traces come from 
the fact that the charmonia are SU(3) singlets. However, one can easily show that 



Thus, there is no electromagnetic contribution to the one pion emission transitions at order 
0(a). Actually, there is a more general relation 



That means, for any transition with the emission of one soft pion between two iso-singlets, 
the contribution from virtual photons exchanged between light quarks vanishes at tree- 
level. 

2) The virtual photon is exchanged inside the heavy quarkonia. In this case, no operator 
containing light mesons without derivative or quark mass can be constructed. This may 
be understood as virtual photons exchanged inside the heavy quarkonia cannot contribute 
to the isospin breaking transitions. 

3) The virtual photon is exchanged between a heavy (anti-)quark and a light (anti-)quark. 
This kind of virtual photon is important in understanding the isospin mass splitting of 
heavy hadrons [63], 164] . In principle, this will give a non- vanishing contribution to isospin 
breaking transitions. For the transitions considered here, however, only one Goldstone 
boson is emitted. Therefore, the operator for the light flavor part should be {Q~) a a- 
However, the trace of Q- vanishes (see Eq. ([22]) ). 

Therefore, there is no e.m. contribution to the isospin breaking heavy quarkonium transitions 
at order 0(a) and thus they can be neglected compared to the quark mass difference terms. 
This conclusion agrees with those of earlier studies in Refs. [331 134] , 

2.3 Tree-level amplitudes 

Before working out the tree-level amplitudes using the Lagrangians given in Section 12.11 
one subtlety needs to be addressed. In Eq. ([3]), the tt° and rj are SU(3) flavor eigen-states. 
However, they are not exactly the same as the physical pion and eta which are mass eigen- 
states. Denoting the physical states by tt° and 77, the ir° — rj mixing is given as 



{Q~)aa = 0, 
(Q+Q-)aa = 0. 



(22) 




(23) 



7r = 7r cos e w o v — fj sin e n o v = 
rj = fj cos e n o n + 7T° sin e^ v = 



7T° - e n o v f) + O (e^oj , 
fj + e^TT + O (el ^j , 



(24) 
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itiAjj k e i (i//)e>(J/il>)q k B du 
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> J/lpT] 
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" Xc07T° 




XcO - 


XclTT° 


-2^/6i7e(x c i) ■ (f-Bdu 


Xcl - 


-> Xcivr 




Xcl - 




3v / 2i7£ i (x / c i)e ii (Xc2)g i 5 dll 


Xc2 ~ 


"> Xc2^° 


-i6^ k e il {x' C 2^ jl (Xc2)q k B du 


K- 


► M° 


-iQ^e^h'^ih^B^ 



Table 2: Tree-level amplitudes for the charmonium transitions with the emission of one pion 
or eta. A factor of y/MiMf, with M^in denoting the mass of the initial (final) charmonium, 
should be multiplied to all the expressions to account for the non-relativistic normalization 
of the charmonium fields. 



where e^o^ is the well-known 7r° — 77 mixing angle, which reads to LO in the chiral expansion 

V3m d -m u 

Mr; = -j — (25) 

' 4 m s — m 

with rh = (m u + md)/2 the average mass of the u and d quarks. Using Dashen's theorem |65| . 
one may express the mixing angle in terms of the masses of the Goldstone bosons at LO in 
CHPT 

1 Ml - ML + M 2 + -M 2 

6,0. = ^ S „,2 = °- 0L 26 

n VS M 2 -M 2 

The mixing of the 7r° or 77 with the rf is not considered since it is of higher order. The reason 
is that the rf is not a Goldstone boson of the spontaneous breaking of SU(3)lxSU(3)r to the 
vector subgroup SU(3)y and its mass as a large scale provides a suppression. 

The tree-level amplitudes for the charmonium transitions with the emission of one pion 
or eta are listed in Table [2^1 where we have defined B du = B^rrid — m u )/F and B s i = 
Bo(m s — m)/F. A factor of yjMiMf, with Mufs denoting the mass of the initial (final) 
charmonium, should be multiplied to all the expressions to account for the non-relativistic 
normalization of the charmonium fields used in the effective Lagrangians. 



3 Decay amplitudes from charmed meson loops 
3.1 Charmed meson loops 

In this section, we list all the possible loops (i.e. the triangle graphs) with the lowest-lying 
pseudoscalar and vector charmed mesons for each transition. There are three charmed mesons 
in each loop. To be specific, we denote the one connecting the initial charmonium and the 

3 We have checked that the ratios among the spin-averaged absolute square of the transition amplitudes for 
the transitions between P-wave charmonia given in Ref. [5S] can be reproduced. 
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(d) (e) (f) 

Figure 3: The decays ip 1 — > J/tpTT°(r]) through triangle charmed-meson loops. Charmonia, 
light mesons, pseudoscalar and vector charmed mesons, are denoted by double, dashed, thin 
and thick solid lines, respectively. 



light meson as Ml, the one connecting two charmonia as M2, and the one connecting the 
final charmonium and the light meson as M3. The meson Mi has a mass m^. For instance, 
in Fig. [3] which shows all loops contributing to the i/j' — > J/i/m°(r]), Ml, M2 and M3 are the 
D, D and D* , respectively in diagram (a). All the loops contributing to each decay are listed 
in Table [3j 

3.2 Leading order effective Lagrangians 

In order to calculate the leading contributions from the charmed meson loops, we need the LO 
effective Lagrangians for the couplings. Because the pion and eta are pseudo-Goldstone bosons 
of the spontaneous chiral symmetry breaking of QCD, their coupling to the charmed mesons 
in the low-energy limit is constrained by chiral symmetry. The effective Lagrangians were 
constructed considering both the heavy quark symmetry and chiral symmetry in Refs. |66} 
[67l [68] (for a review, we refer to Ref. [69]). In the two-component notation of Ref. [56J, the 
charmed mesons are represented by H a = V a ■ a + P a , with V a and P a denoting the vector and 
pseudoscalar charmed mesons, respectively, a are the Pauli matrices, and a is the light flavor 
index. Explicitly, one can write P a {V a ) = (p^°, L> ( * )+ , L>i* )+ ) . The lowest order chiral 
effective Lagrangian for the axial coupling is [56J 

= (HlH b ff-u ba ), (27) 

where the axial current is u = —y/2d(/)/F + O(cj) 3 ). 

The LO Lagrangian for the coupling of the S- or P-wave charmonium fields to the charmed 
and anti-charmed mesons can be constructed considering parity, charge conjugation and 
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v/ - 




[£>,A-D*], 


D,D*,D*l [D*,D*,D], [D*,D,D*], [D*,D*,D*\ 


%p - 


> h c TT° 


[£>, [£>*,£,£>*], 


[D*,D*,D], [D,D*,D% [D*,D*,D*\ 


h c - 


y J/^7T° 


[£>*,£,£>], 


[D,D*,D*\, [D*,D*,D], [D*,D*,D*} 


V'c ~ 


' Xcovr 




[D*,D*,D*\ 


XcO ~ 


->■ ry c 7r° 


[D,D,D% [D*,D*,D], 


[D*,D*,D*\ 


XcO - 




[D,D,D*], [D*,D*,D] 




Xcl - 


-> Xcovr 


[D*,D,D], [D,D*,D*\ 




Xcl - 


-> Xcivr 


[D*,D,D*\ 




Xcl - 


"> Xc2TT° 


[D,D*,D*\ 




X' C 2 - 


"> XclTT 


[D*,D*,D] 




Xc2 ~ 


"> Xc2VT° 


[D*,D*,D*\ 




K- 


► /i c 7T° 


[D*,D,D*\, [D,D*,D* 


, [D* , D* , D], [D*,D*,D*\ 



Table 3: All the loops contributing to each transition. The mesons are listed as [Ml, M2, M3}. 
Flavor labels are dropped for simplicity. 



spin symmetry. In two-component notation, the one for the 5-wave charmonia J/ij) and 
rj c reads [32] 

C^=i^-^H a a-^H a )+R.c, (28) 

^ y _ -* _ 

where A d B = A(dB) — (dA)B, and H a = —V a -a+P a is the field for anti-charmed mesons [57] . 

The Lagrangian for the P-wave charmonia at LO is [57] 

C x = i| (^'Hao'Ha) + H.c. (29) 

These Lagrangians were introduced in Ref. [70] in four-component notation. The values of g\ 
and 52 are twice of those in Ref. [70] • @ The Lagrangians for the coupling of the radial excited 
charmonia to the charmed and anti-charmed mesons have the same form as Eqs. (|28|29p with 
the coupling constants changed to those for the excited states g' 2 and g[. For later use, we 
evaluate the traces in Eqs. ()28|29p . and rewrite the Lagrangians in terms of the meson fields. 
The Lagrangian for the J /if) and rj c is 

= i92^ (yffivi - vffivj - Vffivj) +g 2 e ijk ^ (P a ^ j V a k - vffip a ) 

+i 9 2^Pa^Pa + 92vU ijk vf3 j V* + i 92 vl (v^P a - P Wj) + H.C. (30) 

The Lagrangian for the XcJ and h c reads 

= i9iX ] M {ViVi + ViVi) + V2g lX ] c \ [ViP a + P a Vt) + A=9ixlo (Va • V a + 3P a P a ) 
-gi^hl'V'Vj! + i 9l ht {v:P a - P a V^ + H.c, (31) 

4 In the definition of the Lagrangians of Ref. [70] , each term is doubled for heavy quarkonia with the same 
flavor of quark and anti-quark. Hence the values of the coupling constants there should be half of those in our 
paper. 
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where the trace and symmetry properties Xc2^ = an< ^ Xci^^ = have been used in the 
derivations. 

3.3 Decay amplitudes of the pion or eta emission transitions 

The amplitudes for all the transitions with charged charmed meson loops are listed in Ap- 
pendix [Bj The charmonia are isospin and SU(3) flavor singlets, the pions form an isospin 
triplet and the rj is an element of the SU(3) octet. Therefore, the transitions between two 
charmonia with the emission of one pion or eta break isospin or SU(3) symmetry. The lead- 
ing contributions to the eta transition amplitudes are given by the differences between the 
non-strange and strange charmed meson loops. The decay amplitude for the tp' — > J /ipr] is 

M -> J/^ri) =-j=[M -> W) Q + M (V' -> J/W)± ~ -)• W) J , (32) 

where <p denotes the light pseudoscalar meson, and the expression of Ai (ip' — > J/ip(ft)± is 
given as Eq. (IB.lh in Appendix |Bl The amplitude with lower index can be obtained by 
simply replacing the charged charmed mesons in Eq. (|B.1|) by the neutral ones, and the lower 
index s denotes the charmed-strange meson loops. The leading contributions to the pionic 
transition amplitudes are given by the differences between the neutral and charged charmed 
meson loops, and also from the the ir° — rj mixing through the loops contributing to the eta 
transition. From Eq. (|24p . the physical pion is 

vf = vr + e n o vV + O (4 J . (33) 

Hence, the amplitude for the tp' — > J/tpTr° reads 

M (tp' -> J/ipTT°) = ^p- \M (tp' -> J/4>(f) + M{ip' -> J/tp(p), - 2M (tp' -> J/ip(p) , 
V3 L 

+M {iP' -»■ J/ipcP) Q — M(ij)'-t J/ip(p) ± • (34) 

Note that although the loop amplitudes M. (tp' — > J/ij)<f>)o± s take the same form as those 
in Eq. (|32p . the momentum q is different for the different decays, it is given by the three- 
momentum of the light meson in the final state. For all the other transitions considered in 
this paper, the available phase spaces only allow the emission of one pion, and the amplitudes 
can be obtained in a similar way to Eq. (|34p using the equations given in Appendix [Bl 

As can be seen from the presence of a Pauli matrix between two charmonium fields in 
the Lagrangians Css an d £sp> the transitions between the S- wave charmonia violate the 
heavy quark spin symmetry, so do those between one S and one P-wave charmonia. Were the 
heavy quark spin symmetry exact, the vector and pseudoscalar charmed mesons would have 
the same mass. In the heavy quark spin symmetric world, all the meson- loop amplitudes for 
the transitions between two 5"-wave charmonia, and those between one S and one P-wave 
charmonia would vanish. This is because the contributions from different loops would cancel 
with each other completely as one may easily see from the amplitudes listed in Appendix IB1 by 
putting Mb = Mo* ■ This means the vector and pseudoscalar heavy mesons have to be consid- 
ered simultaneously to keep the structure of the spin symmetry which is used in constructing 
the Lagrangians and relating the coupling constants for different transitions. The transitions 
between the P-wave charmonia respect the spin symmetry, so the resulting amplitudes do not 
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vanish in the given heavy quark spin symmetric world as shown in Appendix IB. 31 Keeping 
spin symmetry structure of the tree-level amplitudes might be a general feature of the heavy 
hadron loops, so that the analysis of the spin partner of heavy hadron molecules in Ref. |71j . 
which is based on the spin symmetry without considering loops, would not be affected. 



3.4 Power counting of the loops 

Before proceeding to numerical calculations, one must analyze the power counting of the loops 
in the NREFT formalism. Any loop diagram in the paper is composed of two vertices for 
the coupling of a charmonium to charmed and anti-charmed mesons — both characterized 
by a vertex structure and a coupling constant, one vertex for the coupling of a pion or eta 
to charmed mesons, and three propagators for charmed mesons. The power counting of a 
given diagram is obtained via estimating each individual ingredient by a typical value. Based 
on this analysis, each diagram can be assigned a definite order n in the given expansion 
parameter, which for the transitions at hand turn out to be the velocity v of the intermediate 
heavy mesons. Once a complete calculation up to order n is performed, the uncertainty 
of the calculation may be estimated as ?/ n+1 ) — we will show below that additional scales 
introduced by the dimensionful coupling constants do not distort this picture. 

Since effective field theories are in general non-renormalizable, they have to be regularized 
and renormalized order by order. A consistent power counting thus has to guarantee that at 
each order where there are divergences there are also appropriate counterterms available. In 
this section we check this for the transitions and the power counting at hand. 

Let us first focus on the propagators, and take the scalar loop integral as an example. The 
scalar loop integral in d dimensions for the triangle graphs under consideration is defined as 

T( W [ ddl i - [JlI 1 

{Q) J (27r) d D J (27r) d (P -m\ + ie) [{P - I) 2 -m 2 + ie] [(I - q) 2 - m§ + ie] ' 1 ' 

where P is the momentum of the initial charmonium and q is the momentum of the light 
meson in the final state. Non-relativistically, in the rest-frame of the initial charmonium, the 
loop can be written as 



dH 1 



Smxmamgy (2vr) 4 fp_ Tl ^ (po_p_ T2 ^ fto_gO_ r3 (|f_g|) 
1 f dH 1 



8m 1 m 2 m 3 7 (2tt)3 ( Ei -T 2 ^y Tl {\l\)) ( Ef -T 2 {\j\)-T z {\l-q\) 



(36) 



where Ti(p) = p 2 /2nrii = rriiV 2 /2, with v being the charmed meson velocity, denotes the 
kinetic energy for the charmed mesons with masses mi, 777-2 an d 777-3, P*i = M% — mi — 777-2 
and Ef = Mt — 777-2 — m 3 — E n denote the energies available for the first (before the pion 
emission) and second (after the pion emission) two-heavy-meson intermediate state. One 
may assign the charmed meson momentum as Mjjv. Here, we will only count the power of 
v since the dimension of the loops can be simply implemented by multiplying proper power 
of Mjj. Thus, the scalar loop scales as 1/(16?t) [v 3 /(v 2 ) 2 ] = 1/(167ttj), since in the last line of 
Eq. ([36]) each of the non-relativistic propagators is counted as 1/tj 2 and the integral measure 
is counted as 7j 3 /(16?t), where it was used that the loops are dominated by the unitarity 
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cut, which produces a factor of ir to be combined with the standard factor l/(47r) 2 from the 
integral measure. This factor is common to all loop contributions as given in Table 01 

As indicated in the Eq. (|27|) for the axial-coupling of the pion or eta to the charmed 
mesons, the corresponding vertex is proportional to the external momentum of the pion or 
eta, denoted by q — this gives one power of q in the expressions for all the loop contributions 
in Table HI Further, we have to account for the isospin or SU(3) violation as well as the 
momentum dependence of the charmonium-charmed meson vertices — the scaling of the 
coupling constants will be discussed below. To account for the corresponding symmetry 
breaking, in each power counting estimate for the loops listed in Table HI we have pulled 
out the meson mass difference, denoted as which is a small energy scale, and divided A 
by a factor which characterizes the intrinsic energy, v 2 , for balance. The vertices are more 
complicated since their scaling behavior depends on the quantum numbers of the charmonia 
in both the initial and final states. We therefore classify the vertices into three groups: 

1) Both the vertices for the initial and final charmonia are in S-wave. This corresponds to 
the transitions between two P-wave charmonia (denoted by PP in Table HJ). In this case, 
the vertices do not give any non-trivial contribution to the power counting, and they scale 
as O(v°q ) = 0(1). 

2) One vertex is in S-wave, and the other one is in P-wave. This corresponds to the transitions 
between one S- wave and one P-wave charmonia (denoted by SP in Table U]). In this case, 
the loop momentum must be contracted with the external momentum q, and hence render 
the scale of the P-wave vertex to be 0(q). In this case a factor needs to be introduced 
to match dimensions. 

3) Both the vertices are in P-wave. This corresponds to the transitions between two S'-wave 
charmonia (denoted by SS in Tabled]). In this case the loops are tensor loops, which 
can be split into two parts as given in Eq. (|A.5j) . The part q l q 3 lQ (q) scales as two 
powers of the eternal momentum, i.e., 0(q 2 ). In the other part, the two momenta in the 
numerator of the loop integrand contract with each other, and as a result, the Kronecker 
delta appears. Here all momenta appearing are internal momenta, which, by assumption, 
scale as v. Thus, for this piece of the integral the two vertex functions together scale as v 2 . 
In the transitions considered in the paper, one always has q < Mpv. Hence, the product 
of the two vertices in this case can be counted as 0{y 2 ). 

The last ingredients to be discussed are the coupling constants g, g\ and gi. g is the 
axial coupling constant for the heavy mesons. It is dimensionless, and should be of order 
unity. Based on the underlying Lagrangians, e.g., the dimension of gi , the coupling of the 
S'-wave charmonia to the open charm ground states, is —3/2. In Ref. [70], using vector meson 
dominance arguments, the authors obtain gi = ^/Mj/^/iMnfj/^), where fj/^ is the decay 
constant of the J/ifi. On the quark level it scales with the J/tp wave function at the origin 

3 II 

that, on dimensional grounds, should be /ja/> ~ m c v c — the quark mass and velocity m c 
and v c were introduced at the end of Section 12.11 Hence, we have 

{m c v c ) 6 ' z 



5 These meson mass differences are, of course, generated by quark mass differences and e.m. effects. For 
the charmed mesons, A is of similar size as the quark mass differences 5. 
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Using m c = 1.5 GeV and v 2 = 0.3 we get gi = 1.9 GeV 3 / 2 , which is close to independent 
model estimates in a range of 2.1. ..2. 9 GeV -3 / 2 for this quantity existing in the literature |85l 



The expression for g\ derived from vector meson dominance is g\ = — 2-y/ M Xc0 /3/ f Xc0 [70]. 
From dimensional analysis, the decay constant of the P-wave charmonium XcO should scale 
with the first derivative of the wave function at the origin. Hence, 

/2 2 
3 — 5/2 ' ( 38 ) 

Using the above scaling, we can work out the power counting of the loops for different 
processes. For the SS transitions, as mentioned at the end of Section [33j the loop amplitudes 
vanish in case of = Mp*, i.e., the loop amplitude violates spin symmetry as the tree level 
amplitude does. So a factor of AqcdA^c should also be considered in the scaling of the loop 
amplitude as well as the estimate for the loop integral itself, l/(167ru), and a factor v 2 , which 
originates from the two decay vertices. Furthermore, the nonrelativistic normalization factor 
MiMf-mimim^ gives a factor 2m^. As explained in the above, a factor of A/(m,DV 2 ) ~ 
A/(m c v 2 ) should also be introduced to account for the isospin breaking. Collecting all factors, 
we get for the scaling of the loop amplitude for the SS transitions 



v 2 Aqcd A 
b Ibivv m c m c v z 



1 1 qA 

47tv 3 m c v 



(39) 



As will be shown later, these rules for power counting of the loops are satisfied by explicit 
non-relativistic calculations. 

Similarly, for the SP and PP transitions, using the scaling of g\ given in Eq. (|38p we have 

KA SP 9 r, 4 1 1 A QCD A 

1 q 2 A 



(40) 



and 



pp g . a 1 A 
^loop ~ F 9mq2m cT6^v— 2 



1 1 qA 



3ttv^ Aqcd v 3 



(41) 



In the last equation, a factor of 1/M|, ~ l/ m c was taken into account to give the correct 
dimension of the scalar loop. This can be done as mentioned below Eq. (|36p . 

The scaling for the loop amplitudes in Eqs. (|39l) . (|40l) and (j4T1) need to be compared to the 
corresponding estimates for the tree-level amplitudes. The relevant estimate is given explicitly 
in Eq. (]18|) for the SS transitions and Eq. (I19p for the PP and SP transitions — the final 
results for all transitions are given in Table HI For charmed mesons, A is of similar size as 5 
as may be checked numerically, i.e. M D + — M D o ~ — m u and M n + — M D + ~ m s — uid. 
At this stage, let us make a remark about the factor containing the charm quark velocity 
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ss 



SP 



PP 



Tree-level 



1 



A, 



QCD 



Loops 



N_ 

m r , 



v 6 \ 2 ( A 



A/" gA 



m r v 



qq_ 

Ml 



> 3 M 



-A 



D 



A, 



QCD 



N qA 



A, 



QCD 



Table 4: Power counting of the tree-level amplitudes and the (leading) loops. Here SS, SP 
and PP represent transitions between two .S'-wave, one S'-wave and one P-wave, and two 
P-wave charmonia, respectively. The parameter 5 denotes the quark mass differences, and A 
the charmed meson mass differences. They are the strength parameters for isospin or SU(3) 
symmetry violation, v is the heavy-meson velocity in the intermediate loops, q the momentum 
of the outgoing pseudoscalar meson, and Md the mass of the heavy mesons in the loop. For 
the origin of the individual factors, see the text. TV = 1/(4tt^), 1/(2^/3^) and l/(37rof) 
for the SS, SP and PP transitions, respectively, where v c denotes the charm quark velocity 
inside the charmonia. 



scaling of the coupling constants, i.e. N = 1/(4^^), l/(2v / 3vrf^) and 1/(3ttv^) for the SS, 
SP and PP transitions, respectively. Taking v c = V0.3, we get N = 0.5, 1 and 2 for the 
SS, SP and PP transitions, respectively. All these numbers are of order unity. Although J\f 
scales differently in powers of v c for different processes, the numerical values are not very 
different. This is because the charm quark velocity v c in such a charmonium is not so small. 
Hence the numerical values of the presumably well-separated scales may be similar in practice. 
Sometimes, the order of the scales is even reversed. For instance, purely from the scaling, 

3/2 5/2 

fj/ip ~ fncVc should be larger than f Xc0 ~ m c v c ■ However, using the experimental value 
r(J ' /ip — > e + e~) = 5.55 ± 0.14 keV [Uj, and the relation between decay constant and the 
leptonic width of the J /tp 

r(W- + = ^_4/V (42) 

fj/tp = 416 MeV, while the numerical result from QCD sum rules gives f Xc0 = 510 ± 
40 MeV |72j . Furthermore, in view that there must be unknown numerical factors in the 
coupling constants g\, 52 and also the the tree-level ones, and the v c scaling of tree- level 
couplings is not known yet as mentioned before Eq. (fl~6|) . we will neglect the subtlety caused 
by M in the power counting, and just take it to be unity for all the transitions (numerical 
calculations in the next section will not be affected). Based on the underlying power count- 
ing, we therefore predict loops to be enhanced by a power of 1/v for SS and l/v s for PP 
transitions, while for SP transitions in many cases loops appear to be suppressed — see 
detailed discussion below. One should also keep in mind that the meson velocity v in some 
of the processes to be discussed is as large as 0.5 so that 1/v for the SS transitions might 
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be in practice not a large enhancement. But it should be sufficient to say that the the loops 
are important compared to the tree-level contributions, and should not be neglected in any 
realistic analysis. 

As mentioned above, a consistent power counting needs to ensure that all appearing di- 
vergences can be absorbed into appropriate counterterms. This is of importance here since 
the leading counterterms are supplied by the tree-level amplitudes discussed above, which in 
some cases appear in higher orders than the leading loops. In such a case, they might get 
renormalized by absorbing the divergence of the loops. Although the scalar loop, defined 
in Eq. (I35j) is convergent, for some transitions there are momentum factors at the vertices, 
resulting in divergent integrals. In Appendix [B] the contributions to all transitions are ex- 
pressed in terms of a few fundamental integrals defined in Appendix El Note that, in order to 
account for the non-relativistic normalization of the charmonium and charmed meson fields, 
one needs to multiply each amplitude by a factor of y / M^M fm\m2m^, where M^n is the 
mass of the initial (final) charmonium, and nrii (i = 1,2,3) are the masses of the charmed 
mesons in the loops. The loop functions I^(q), Iq (q) and l[ 2 \q) are constructed from the 
basic scalar loop functions I(q) defined in Eq. (l35j) 



= ' I wfd ■ + - ' I • (43) 

and B(c), which in d space-time dimensions is 

B(c) = AttN / " . , - , (44) 

V ' J (2vr)^i/2 +c _ ie ' v > 

where N = ^i2^23/(16vrmim2m3) with fiij = rrnmj / '(m, + rrij) for the reduced masses. For 
later use, we define the following quantities 



M23 \ _ 2 o > / > , ^23 2 
) q , C= 2^i2&12, C = 2/i 2 3»23 H q , 

m 3 J m 3 



(45) 



where 612 = mi + m 2 — Mj, and 623 = ^2 + "*3 + q° — M^. 

In our full calculation all integrals are evaluated using dimensional regularization. In 
this section, however, in order to make all divergences explicit, we investigate the divergence 
structure by introducing a sharp momentum cut-off in three dimensions. With this one finds 

The non-analytic part is finite. It is determined by unitarity, and hence does not depend 
on the choice of regularization method. In the integral 1^ only the linear combination 
B(c' — a) — B(c) appears — c.f. Eq. (|A.6|) — and thus the divergence of B(c) cancels. The 
UV divergent part of q 2 l'^'\{q) is 

flr 2 J< 2) (g) uv = r 2 Ii (a) (ff) w = -A. (47) 

7T 

For the loop amplitudes given in Appendix [Bl only the transitions between two S'-wave 
charmonia, i.e. ip' — > J j 'ipir (77) , are divergent. 
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In order to estimate the finite parts of the pertinent integrals we may use a -C c ~ d ~ 
\J2[ib ~ Mdv and expand the expressions in a series around a = 0. With this we find 

5(0' - a) = 5(c) + 0(a) = -AVc + 0(a), 
/(g) « iV_L + 0(a), 

/«(,)« A^J- + 0(a), 

^2 7 (2) (9) finite ^ .jy^ + 0(a), 

^2 7 (2) (g) finite = 0(fl) (4g) 

The exact expressions are given in Appendix 1X1 

As an example we now focus on the analysis of diagram (b) in Fig. [3] to the amplitude of 
the if/ — > J/V'tt — the discussion is easily generalized to the other diagrams. The amplitude 
reads 

MM J/V7r°) (b) = N (h) [q 2 l{ 2 \q,D* ,D ,D°)-q 2 l[ 2 \q,D* ± ,D ± ,D ± )] , (49) 
with Nr b \ = —4(g/F)g2g' 2 ( l: > k q t e^(Tp')e k (J/Tp). The contribution of the finite part of the loop 

(2) 

function I\ (g) to diagram (b) behaves as 

M (V' -»• J/^°)^f e ~ -%) (iV n V2^X - A^c V^lX) , (50) 

where the lower index n means neutral, and c charged. Denoting the mass difference between 
the charged and neutral charmed mesons by A@ we have // c = \x n + A/2 and 6 C = 6 n + 2A. 
Thus, we have 

Mtf J/V)®* ° ~ jV (b) ATA 2/i - + 6 " /2 + 0(A 2 ) 

~ JV(b)^~, (51) 

which is consistent with the power counting analysis for iSS" transitions given above. Here, 
we neglect the difference between N n and N c since it is of higher order. On the other hand, 
the UV divergence of diagram (b) is 

MW J/^tt )^ = N (h) ± (N n - N c ) 

A A 

~ N (h) N . 52 

Therefore, for diagram (b), the finite part is of order (Av _1 ) , while the UV divergence is of 
order (At>°) . Hence the UV divergence is one order higher in the expansion of v. Recalling 
the tree-level contribution to the ifi' — > J/tf^ir starts from the same order as (Av°) - 



We neglect the difference between the mass difference for the pseudoscalar mesons and that for the vector 
ones, which is of higher order in heavy quark expansion. Empirically, one finds Ap = m D ± — m D o = 
4.77 ± 0.10 MeV, and A v = m D ,± - m D ,o = 3.29 ± 0.21 MeV. Their difference is about 30% of A P , which 
can be understood as C(Aqcd/tIc) effects. 
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see the column for SS transitions in Tabled] — such a divergence can be renormalized by a 
counterterm in the Lagrangian for the tree-level contribution. 

As a result of this analysis we summarize, comparing the loop contributions with the 
tree-level decay amplitudes given in the last section and in Table |3J where we assume the 
same scale for the light quark mass and the heavy meson mass differences, 5 ~ A |27| . that 
the loop contributions for the SS transitions are enhanced by a factor of 1/v ~ 2, and for the 
PP transitions even by a factor of 1/v 3 ~ 10. The situation for the SP transitions should 
be analyzed case by case since an enhancement factor 1/v 3 competes with a suppression 
factor q 2 /Mp. For the SP transitions with small phase space, the external momentum q 
might be small enough to make q 2 /(v 3 M^) much smaller than 1, which is satisfied for the 
decays tp' — > h c ir° and r( c —> XcO 7 ? [48J. In this case, the decay is dominated by the tree-level 
contributions. However, for the decays with external momentum q > Mpv 3 / 2 , the factor 
q 2 /(v 3 AIp) > 1, so is no more a suppression, and hence the tree-level contributions are at 
least as important as the loop contributions. In summary, our power counting was shown to 
be consistent with the divergence structure of the pertinent integrals. 

4 Results for the decay widths 

In this section, we give the results from explicit calculations of all the mentioned transitions 
with emphasis on the contributions from charmed meson loops. In the numerical evaluations 
we use the following values for the meson masses [44] 



M w o 


= 134.98 MeV, 


M v = 


547.85 MeV, 






M D+ 


= 1869.60 MeV, 


M D o = 


: 1864.83 MeV, 


M Ds = 


1968.47 MeV, 


M D ,+ 


= 2010.25 MeV, 


M D * -- 


= 2006.96 MeV, 


M D * = 


= 2112.3 MeV, 




= 3096.92 MeV, 


M Vc = 


2980.3 MeV, 


My = 


3686.09 MeV, 




= 3637 MeV, 


M hc = 


3525.42 MeV, 


M Xc0 = 


3414.75 MeV, 




= 3510.66 MeV, 


M Xc2 = 


= 3556.20 MeV, 


M x> 2 = 


= 3929 MeV. 



(53) 



Here we have identified the x' C 2 with the Z(3930) as done by the Particle Data Group, see 
e.g. Refs. [73l EH El [50] . The Z(3930) was observed in the DD mass distribution in photon- 
photon collisions by the Belle Collaboration [73] with a mass of 3929 ± 5 ± 2 MeV and width 
of 29 ± 10 ± 2 MeV. The observed angular distribution suggests that its quantum numbers 
are J PC = 2 ++ '. The widths of the x' C 2 calculated in quark models are consistent with the 
observed values for the Z(3930) [751 1511] . 

There is no unambiguous candidate for either of the other three excited P-wave charmonia 
x' c0 , x'd an d h' c . Very recently, the Belle Collaboration observed an enhancement, called the 
V(3915), in the J/tpco mass distribution in photon-photon collisions [75J. The mass and width 
of the X(3915) were reported to be 3915±3±2 and 17±10±3, respectively [75]. In Ref. [76], 
the authors suggest the X(3915) as the x' c o- However, this assignment may be criticized from 
the following points: First, since the x' c o 1S above the DD threshold, and it couples to DD 
in an S- wave, one would expect it to have a larger width than the 17 MeV of the X(3915). 
Second, were the Z(3930) the x' C 2> the mass difference -^^(3939) — M X (3gi5) = 14 =t 6 MeV 
is too small for the hyperfine splitting, since the mass splitting between the ground state 
Xc2 and XcO is one order of magnitude larger, M Xc2 — M Xc0 = 141 MeV [UJ. In addition, 
one expects the hyperfine splittings for the bottomonia are smaller than the corresponding 
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ones for the charmonia, as can be checked from all the measured cases. The mass splitting 
between the excited P-wave bottomonia x'b2 ano - Xbo * s 36.2 ± 0.8 MeV [Uj. It is larger 
than M^(393o) — Mjf(39i5) 3 and therefore does not support the assignment of the A"(3915) and 
Z(3930) as the x' c o an d x' C 2 simultaneously. Based on the above arguments, we shall let the 
mass of the x' c o run i n a range from 3800 MeV to 3930 MeV which covers the predicted values 
from quark models [TTJ EH] • 

In the observed spectrum of the charmonia, the only candidate of the x'& with well- 
established quantum numbers is the X(3872) discovered by the Belle Collaboration [78] . 
However, due to the proximity to the DD* threshold, the interpretation of the X(3872) as a 
molecular state [79l[HD] or virtual state [81] (for an update of the latter analysis see Ref. [82] ) 
is very intriguing (for reviews, see, e.g. Refs. El [8] ) . Therefore, although a molecular 
interpretation of the X(3872) was questioned in Ref. [83] (see also Ref. [84] for a critical 
re-evaluation), we shall not identify the X(3872) as the x!cii an d the mass °f the x'& wm a ^ so 
be allowed to vary. To be specific, a range from 3.83 GeV to 3.93 GeV will be chosen which 
covers the the predicted value from some quark models, e.g. Ref. [50J, and the mass of the 
V(3872). 

4.1 Transitions between the S-wave charmonia 

The decays studied in this class are ip' — > J/ipir° and ip' — > J/iprj. As shown in Sec. 13.41 (see 
especially Table EJ), for these transitions the charmed meson loops are enhanced by 1/v. The 
velocity can be estimated as v ~ <J '(2M^ — M^)/M^ ~ 0.53 with being the averaged 

charmed-meson mass, and = {Mj/^, + Mfi)/2. Hence, the LO result for the width 
is provided by the loops [32], although the relatively large expansion parameter leads to a 
sizable uncertainty. 

Taking into account only the loop contributions, we give the numerical results for the 
ip 1 — > J/^7r° and ip 1 — > J/tpij. To account for the non-relativistic normalization, one needs to 
multiply the amplitudes by a proper factor. 

One way to do this is to multiply the amplitude for each loop by the factor ^fM^MjJ^ x 
mim2m3, with mi being the masses of the charmed mesons in the loop. In this way, the 
factor may differ for different loops. From the decay D* — > Dtt the coupling g, defined in 
Eq. (|27p . can be fixed to g = 0.6 using F = 92.4 MeV. We thus get for the decay widths 

T(V/ -> J/V>tt°) = (0.048 ±0.025)^2 2 keV, 
r(V' J/ih) = (0.43 ± 0.23)g 2 2 g' 2 2 keV, (54) 

where an uncertainty of 53%, which is the value of v, has been taken into account. Here and 
in the following values for g 2 and g' 2 are given in unit of GeV -3 / 2 . The resulting ratio, which 
is parameter-free, reads 

i^o^ = 0.11 ±0.06, (55) 

with an uncertainty of 53%. Comparing with the measured values listed in Tabled! it is within 
two sigma of the CLEO and PDG-fit data, and even overlaps within uncertainties with the 
result given by the BES Collaboration. A more conservative estimate of the uncertainty of 
the ratio may be given by assuming the uncertainties in Eq. (|54|) are uncorrelated, which 
would give R^o/^ in a much larger range from 0.03 to 0.36. It is consistent with the data in 
Tabled) 
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Within the effective field theory we can not predict the absolute rates for the decays, for 
the couplings are unknown. However, we may use the results given above to extract some 
averaged coupling constants that may then be compared to the corresponding values in the 
literature. Using the experimental data for two decay widths T(ip' — > J/ip7r°) = 0.40±0.03 keV 
and T(ip' — > J/iprj) = 10.0 ± 0.4 keV |H], the value of G = \fg2g\ can be deduced. Since 
with the more conservative uncertainty estimate a ratio consistent with the data is obtained, 
a combined fit to both widths is possible. The widths for both decays agree with the data 
within uncertainties giving G » 2.0 GeV~ 3/2 . We can define dimensionless coupling constants 
g , (*) (*) = #2 [MjuhM '{*)M ' (*) and similar quantities for the if)'. Assuming the coupling 

constants are SU(6) symmetric, i.e. g^DD = 9i y Mj^m 2 with fh being the average mass of 

the charmed meson spin-flavor SU(6) multiplet, we obtain G^dd = ^9ip D D9ip' D D ~ 7.3. 

Another way to account for the non-relativistic normalization is to use an overall normal- 
ization factor, i.e., UM^i Mj /^fhr . In this way, the results are collected as follows 

r(tp' -> J/ipir ) = (0.098 ± 0.052) 5 | ff2 2 keV, (56) 
r(^' -> J/^) = (0-84 ± 0.45)^^ 2 keV, (57) 

and their ratio is almost the same as that given in Eq. (|55[) . The combined coupling constant is 
extracted from a combined fit as G ~ 1.7 GeV~ 3 / 2 and the dimensionless one is G^od ~ 6.0. 

Note that because the ip' and J ftp are below the open charm threshold, their coupling 
to the charmed mesons cannot be extracted directly using the decay widths. Nevertheless, 
there are several theoretical estimates. Using vector meson dominance, the value for the 
J/tpDD coupling was estimated to be 7.7 in Ref. [85] and 8.0 ± 0.5 in Ref. [86]. Using QCD 
sum rules, Ref. [87] obtained 8.2 ± 1.3, and using an SU(4) chiral model, it was estimated 
to be 4.93 in Ref. [88J. Assuming g^DD = 9ip'DD, the J/tpDD coupling extracted from 
the ip' —7- J/ipTT°(r]) considering only the charmed meson loops is of similar size as these 
phenomenological estimates. 

Since the results for the widths considering these two ways to account for the non- 
relativistic normalization are consistent within uncertainties, in the following, we shall only 
use the former one, i.e. using the physical masses of the intermediate mesons in the factor 

y/MiMf 771-17712777,3. 

4.2 Transitions between the S- and P-wave charmonia 

In these transitions, the situation is different. As analyzed in Section 13.4} the loops do not 
necessarily dominate the transitions. Especially for those decays which have a very small 
phase space, the contribution from the charmed-meson loops is highly suppressed [48] , For 
these decays the momentum of the emitted pion is so small, that q 2 /(v' i M^ ) ) <C 1. Two 
examples are given by ip 1 — > h c 7r°, with q = 86 MeV, and 7]' c — > XcO 71 " ; with q = 171 MeV. 
Considering that the velocity v is about J|2M^ — (Afy + M/ lc )/2]/M^ ~ 0.4, in case of the 
former transition the dimensionless factor is 

1 q 2 

^^V« 0.03. (58) 
v 6 m D 

Thus, if all couplings are natural, the heavy-meson loop contributions are estimated to only 
give a few percent correction to the tree-level contribution. In case only loops are considered 
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we get for the width of the tp 1 — > h c 7r° the very small prediction, 

r(V' -> h c 7T°) loop = 2.1 x W- 7 gig' 2 2 keV, (59) 

where the value of g\ is given in unit of GeV -1 / 2 . The decay width has been measured by the 
BES-III Collaboration very recently [89], and the absolute value of the branching ratio for 
the ip' h c ir° is B(ip' -> ir°h c ) = (8.4± 1.3 ± 1.0) x 10" 4 . Using the PDG value for the width 
of the V', IW = 304 ± 9 keV [43], the measured width is T(i// /i c vr°) = 0.26 ± 0.05 keV. 
Taking into account that g' 2 is about 2 GeV -3 / 2 as extracted in Section \A. II and Ref. [32], and 
<7i is about —4 GeV 1 / 2 from an estimate using the vector meson dominance [70], the numerical 
value T(ip' — > /i c 7r°)ioop ~ 10~ 5 keV is orders of magnitude smaller than the measured value. 
This confirms our very rough estimate for the loop contributions in the above. 

The same situation happens to the rf c — > XcO 71 " - For this decay, the momentum of the 
pion is 171 MeV, and the suppression factor is 



;?2 



v 3 m 2 D 



0.1 . (60) 



Thus, here, for natural couplings, the loops are expected to give a correction of the order of 
10% to the tree-level amplitude. If again only loops are considered, the width is again quite 
small 

r(t/ c ->• Xco Aoop = 1-0 x 10- 5 5l 2 ^ 2 keV. (61) 

The pion momentum is 382 MeV for the transitions h c — > J/ipir , and 387 MeV for the 
XcO ^c 71 " - The dimensionless suppression factor is about 0.3, still a small number although 
much larger than the previous two cases. Especially through the interference with the tree- 
level amplitude meson loops may give a significant contribution for the mentioned decays and 
should not be neglected in any quantitative analysis. The widths, considering only the loop 
contributions, are given by 

T(h c -> J/Vvr°)i oop = 1.9 x 10- 4 g 2 ig 2 2 keV, 
r(x c o -> r? c ^°)ioo P = 3.3 x W- 4 gfg 2 2 keV. (62) 

For the transitions x' c o ~ * Vc^ an d h' c — > ip'ir , the situation is more complicated. Al- 
though the pion momentum is small for these decays, the dimensionless factor {q^/Mn) 2 /v 3 
does not necessarily gives a large suppression. For these two transitions, the masses of both 
the initial and final charmonia are close to the thresholds of the charmed mesons involved in 



the loops. As a result, the velocity as approximated by y \b\/Mi), where \b\ = (\bi 2 \ + \b 2 ^\)/2, 
is small. As discussed at the beginning of this section, the x' c o an d h' c have not been unam- 
biguously identified so far. To illustrate the point, one may use values for their masses from 
quark model calculations. For the x' c o ~^ Vc^ ' considering the loop [D, D, D*], the velocity 
is about 0.30 if we use = 3842 MeV predicted in Ref. [50], and (g,/M D ) 2 /u 3 « 0.23 is 



still a suppression factor. If we use My = 3916 MeV predicted in the quark model [77], the 
velocity is about 0.33, and (q^/Mj^) 2 /v 3 « 0.43. The loops are not highly suppressed, and 
they may have a contribution comparable to or at least not much smaller than the tree-level 
one. Taking into account only the loops, one gets 

r(x' c0 -> ?7cAoop = 0.002(0.010)^^ keV, (63) 
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Figure 4: The widths of the decays x' c o ~~ ^ ^c 71 " an d h' c — >• ^'ir° considering only meson loops. 



where the numbers outside and inside the parentheses are obtained using M x > o predicted in 
Ref. [50] and [77], respectively — the difference reflects the difference in the factor (q w /Mu) 2 /v 3 
which enters squared in the expression for the width. In Fig. 0] (a), the width of the decay 
x'co ~ * Vc 71 ® as a function of My is shown. For the decay h' c — > i/j'tt , considering the loop 
[D*,D,D], the velocity is about 0.14 using My c = 3908 MeV given in Ref. [50]. As a re- 
sult, the factor (q n /M D ) 2 /v 3 « 2.9 even gives an enhancement. Then the width in this case 
is induced mainly by the charmed meson loops. The result for a varying My is shown in 
Fig. [4] (b). In the figure, the thresholds of the DD and D*D* are approached at the lower 
and higher end, respectively. Consequently, the curve for the width shows an increasing ten- 
dency at both ends. Because the width of the x' c o ~^ ls tree-level dominated, and that 
of the h' c — > ip'ir is dominated by loops, which give an enhancement, we predict their ratio 
to be significantly smaller than the one derived from the assumption that both are tree-level 
dominated 

r(x' c0 ^y'y) tvce = M v ,M h , 



25 



0.30 

> 



a 0.25 ; 




~ 0.15; 

^ ; 
>< 0.10 1 

T o : 
~J 0.05 ; 

u 

0.00 r 

3.80 3.82 3.84 3.86 3.88 3.90 3.92 

M(Xa) (GeV) 
Figure 5: The width of the decay x' c o ~^ Xci^°- 



where q\ and q2 are the pion momenta in the rest frame of the initial state for the decays 
Xco ~~ * Vc^ an d h'c ~~ ^ ip'ft ) respectively. Once the masses of the x' c o an d h' c are measured, 
our statement on the ratio can be tested. 



4.3 Transitions between two P-wave charmonia 

As has been shown in Section [3.44 the charmed meson loops scale as qA/v 3 . Compared to 
the tree-level amplitudes, the loops are enhanced by a factor of 1/w 3 , see Table [2J 

Therefore, it is reasonable to neglect the contributions from the tree-level diagrams for the 
transitions between two P-wave charmonia. There are in total seven transition processes from 
the first radial excited P-wave charmonia to the ground state ones. All of the amplitudes are 
proportional to the product of the same coupling constants ggig[/F. Therefore, the ratios 
among these decay widths can be predicted without any free parameter. 

Instead of taking some value of the unknown coupling constants g\ and g[ from model 
estimates, we choose to show the coupling-constant dependent width for the x' c o ~> Xci^° in 
Fig. The curve is obtained by setting the coupling constants g\ and g[ to 1 GeV~ 1//2 , and 
the width can be obtained by multiplying the value plotted in the figure by gfg'i- 

The mass of the x' c i is a ^ so allowed to vary, and the width for the x' c i ~^ XcO 71 " is 
shown in Fig. [6] (a). In Fig. [6] (b), the parameter- free ratio of the decay widths T(x' c i — > 
Xeo7r°)ioop/r(Xco Xci7r°)ioop is shown. In the plots, a double-cusp structure is prominent. 
These cusps correspond to the thresholds of the neutral and charged DD* mesons. When the 
thresholds are approaching, the velocities of the intermediate charmed mesons decrease, and 
then the enhancement characterized by 1/v 3 becomes larger. A similar double-cusp structure 
was found in a study of ao(980) — /o(980) mixing in the decays of charmonia [231 [22 EH], 
where the cusps are due to the thresholds of the charged and neutral kaons [91] . For the 
decays of the x'ci ~~ ^ Xcj^ , the ratios defined as 

ploop _ r(Xcl ~^ Xc07r°)loop ploop ^(Xel — > Xcl^loop 

-"-1012 — fTTi H Z tt0\ ' ^1112 — T^TTj ~ Z tt0\ ' 
1 vXcl ->• Xc2^ )loop 1 Ucl ~^ Xc27T U jloop 
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Figure 6: (a) The width of the decay x'ci ~^ XcO 71 " - (b) The parameter-free ratio between 
the two decay widths T(x' c i ->■ Xco7r )ioop/r(x' c o -> Xci7r )i OO p, here = 3842 MeV as 

predicted in Ref. [50J is used. 



are plotted in Fig. [7] as the solid and dashed lines, respectively, in a parameter- free way. One 
may also construct similar ratios taking into account only the tree-level contributions, which 
are 



p i iv. __ ^InO M Xc0 trce _ 3g^ x M Xcl 

-"1012 — T~3~Y7 ' -"1112 — T - 3~ T7 > V DD J 

5^2 M Xc2 5^2 M Xc2 



where q n j is the pion momentum for the decays x'ci ~^ Xcj^°- The dotted and dot-dashed 
curves in Fig. [7J represent these ratios, respectively. Comparing the curves considering only the 
meson loops with the tree- level ones, the difference is not tremendous. Since the uncertainty 
is large in our loop calculations, c.f. Eq. (|55|) and the discussion below, it may be difficult to 
draw a definite conclusion from the comparison without a refined uncertainty analysis. 

The ratios for the tree-level contributions agree with those derived in Ref. [58J using chiral 
Lagrangians, but differ from those given for the cc option of the X(3872) in Ref. [92] using 
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Figure 7: The parameter-free ratios among the decay widths of the x' c i ~^ Xcj^ considering 
only the meson loops. The corresponding tree-level predictions are also shown. 



the QCDME. However, as shown above, even if the X(3872) is a cc state, its decays into 
the XcJ 71 " are mainly given by the intermediate charmed meson loops, i.e. by non-mulitpole 
contributions. 




M(Xci) (GeV) 



Figure 8: Ratio of the decay widths x' C 2 ~~ ^ Xd^° and x' C 2 ~~ ^ Xc2vr° to F(x' c i ~^ Xc2^ )- 



The widths of the x' C 2 decays into the Xci^° and Xc2^° considering only the meson loops 

are 

r(x' c2 -> Xcl Aoop = (°- 08 ± O.O3)$?0? keV, 

r(x' c2 X c2 ^°)ioop = (0.10 ± 0.04)^f keV, (67) 
where the uncertainties are 33% and 36%, which are the velocities of the intermediate charmed 
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Figure 9: The width of the decay h' c — > h c TT° considering only the meson loops. 



mesons for these two processes, respectively. The ratios of decay widths defined as 



R 



loop 
2112 



r (x' c2 -> Xcivr )i oop loop _ T(x' c2 -> Xc27r°)ioo P 
rC, 



T(Xci ->■ Xc27r°)ioop : 



2212 



r (Xci -> Xc27r°)ioop : 



for the loop contributions only are shown as the solid and dashed lines in Fig. [SJ The same 
ratios considering only the tree-level contributions 



3«3 My , My' 



l 2112 



5q 3 2 M x , c2 M Xc2 



9a? 

ntree ^Hl X c i 
JXn 



•-2212 



5 ?2 3 ^X^ 2 



(69) 



where % is the pion momentum for the corresponding decay, are given by the dotted and dot- 
dashed lines in Fig. [HJ respectively. The difference is obvious, and may be tested by future 
experiments. The width of the decay h' c — > h c 7r° considering only the meson loops is shown 
in Fig. [9j The double-cusp structure is again due to the coupling to the intermediate DD* 
mesons. 



5 Comparison with the effective Lagrangian approach 

As emphasized before, the power counting of NREFT provides a control of the theoretical 
uncertainties. However, for transitions where the mass difference between the initial and 
final charmonium exceeds ~ 800 MeV, i.e. 0(A X ), the NREFT is not applicable any more. 
The ELA, on the other hand, which deals with the non-local effects by introducing empirical 
form factors to cut off the divergences, can be applied to broader kinematic regions. But 
the disadvantage is that, due to a lack of knowledge about the behavior of the counterterms, 
model-dependence will be present in association with the cut-off energies and different forms 
for the form factors. By comparing these two approaches with each other, we expect that the 
mo del- dependent aspects of the ELA can be identified and reduced. In this way, one may 
have more confidence in the ELA when it is applied in kinematical regions where the NREFT 
is not applicable. 
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In this section, we first present the general formulae for the ELA. The calculation results 
will be summarized and compared with the NREFT for those transitions studied in the 
previous section. The detailed formulation is given in Appendix O 

The loop transition amplitudes for the transitions in Fig. [2] can be expressed in a general 
form in the ELA as follows: 

M f*= fp* e — n^^ 2 ) (™) 

J (27r) 4 ^ a x a 2 az 

D* pol. 4 

where V{ (i = 1, 2, 3) are the vertex functions; at = qf — mf (i = 1, 2, 3) are the denominators 
of the intermediate meson propagators. We adopt the form factor, fT- Tiimi, qf), which is a 
product of monopole form factors for each internal meson, i.e. 

Y[Ti(mi,qi) =^1(7711, qx)T 2 {m 2 ,ql)T^{m^ql) , (71) 

i 

with 

Fi(T«,tf)= (^4), (72) 

where A* = rrn + oAqcd [93] and the QCD energy scale is Aqcd = 220 MeV. The numerator 
is introduced in order not to modify the expression at the on-shell point. This form factor is 
supposed to parameterize the non-local effects of the vertex functions and to remove the loop 
divergence in the integrals. In this approach the local couplings for a charmonium to charmed 
mesons, or a light meson to charmed mesons are the same as used in the NREFT, while the 
form factor parameter a will be determined by comparison to experimental information. 
Thus, it is assumed here that all — or at least the dominant part — of the short-range 
physics related to meson loops can be parameterized in the form of Eq. (|72p with a reaction 
independent parameter a. 

Although used widely and very convenient for the actual evaluation of the pertinent 
integrals, the form factor parameterization given in Eq. ()7ip also has its drawbacks. Es- 
pecially, unphysical thresholds, which are located at Aj + rrij = m; + rrij + a Aqcd and 
Aj + Aj = rrii + rrij + 2a Aqcd , are introduced into the integrals. Thus, for sufficiently heavy 
decaying charmonia and small values of a, these singularities are located nearby or even inside 
the physical regime and additionally introduce unphysical contributions. Clearly, those are 
not part of the NREFT. We therefore expect, and indeed observe, that there exist significant 
quantitative differences between these two approaches in the decays of heavy charmonia. It 
suggests that a different form factor parameterization or larger values of a should be consid- 
ered. On the other hand, for those cases, where the form factor singularities do not contribute, 
we expect interesting insights from the comparison of the ELA and NREFT results. 

5.1 Transitions between the 5- wave charmonia 

For transitions between the S- wave charmonia, i.e. ip' — > J/ip7r°(rj), we plot the a-dependence 
of the partial decay widths of ip 1 — > J/i/jTr°(r}) in Fig. \W\ (a) as shown by the solid and dashed 
lines, respectively. The ir° — r] mixing has been taken into account. Using a in a range 
of a ~ 1...2, the measured ratio R^o/^ by both the CLEO and BES Collaborations can be 
reproduced, and the central value of the PDG fit as listed in Table [Q may be obtained with 
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Table 5: Summary of the qualitative and quantitative features in the comparison between the 
ELA and NREFT results for various decay transitions. For PP transitions, we take the width 
ratios in order to compare the ELA with the NREFT results. The appearance of the crosses 
can be understood, as explained in the text. 



a = 1.64, 

r^'^J/^vr ) = Vmiglgi keV, (73) 
r(V/^J/Vr?) = 0.77 52 V keV, (74) 

where the vertex couplings are taken the same as those in the NREFT. These results are 
consistent with the NREFT. 

In Fig. [10] (b), the ratio R n o/ v calculated in the ELA is given in terms of a. The pre- 
dominant feature is that the a-dependence of the ratios is quite stable. It is because the 
loop amplitudes for the transitions play a dominant role in the transition and have the same 
divergence structure which is regularized by the form factors. With a = 1.64, the resulting 
ratio P\oJ v is consistent with the result from the NREFT |32| . 

5.2 Transitions between the S- and P-wave charmonia 

Recall that the power counting suggests suppressions on most of the transitions between the 
S and P-wave charmonia via intermediate meson loops, such as, tp' — > h c ir , while the loops 
in the transition h' c — > i/j'tt° are enhanced. The dominance of tree-level contribution means 
that the physics is described by contact interactions - sometimes loosely called short-distance 
physics. On the contrary, in the ELA, since form factors are adopted, part of the short-range 
physics is already included in the loops. Therefore, one may expect large differences between 
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the NREFT and the ELA in this sector, and the ELA results would typically be larger than 
those in the NREFT, especially for the tree- level dominant transitions. 

With the form factor parameter determined in the previous subsection and the same 
couplings for the vertices as in the NREFT, we now examine the predictions from the ELA 
for those transitions. Instead of going to all the channels, we shall concentrate on two sets of 
decays, i) ifi' — > h c ir and h c — > J/i/m°, and ii) rf c — > XcO 71 " and XcO —> Vc^°, and summarize 
the others in Table [5j 

The two transitions in each set involve similar coupling vertices, but different kinemat- 
ics. Their ratios will again eliminate the uncertainties arising from the unknown coupling 
constants. In Fig. [TTj the decay widths with the charmonium-L'-meson couplings nor- 
malized to unity are plotted for these two pairs of decay channels, i.e. ip' — > h c ir° (solid 
line), h c — > J/ifjir® (dashed line), rj c — > XcO 71 " (dot-dashed line) and XcO — > Vc^° (dotted 
line) in a range of a = 1...2. Note that the BES-III measurement gives a central value 
pexp^/ _^ fi c ir®) = 0.61 keV |89j . Within a = 1...2 which is found in the transitions 
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Figure 11: a-dependence of the decay widths of ip' — > h c 7r° (solid line), h c — > J/t(m (dashed 
line), XcO — > Vc^° (dotted line), and rj c — > XcO 71 " (dot-dashed line). The coupling is G 2 = g\g'\ 
for tp' — > h c ir° and rj' c — > Xc0^°, while G 2 = g\g\ for h c — > J/tpir and XcO Vc^°- 



ip' — > J/ipn° (77) , the partial decay width of — > h c -K° from the meson loops is much smaller 
than the experimental data [89|, which confirms the suppression of the meson loops found 
in the NREFT, although these two approaches give quite different values for the meson loop 
magnitudes. This comparison is listed in Table [5] as qualitative agreement. 

Such a suppression also occurs to h c — > J/i/jtt in the ELA, and exhibits some peculiar 
features in comparison with the NREFT expectation. Since the decay h c — > J/ipn has a 
relatively large phase space, the power counting suppression in the NREFT is not as much 
as for tp' — > h c ir° as pointed out in the previous Section. In the ELA, however, as shown in 
Fig. [TTT these two decay widths are of the same order in the range of a = 1...2. Nevertheless, 
the suppression of the charmed meson loops as analyzed in the NREFT in both decays is 
confirmed in the ELA. The partial decay width of the transition h c — > J/ipir° considering 
only loops is even quantitatively consistent with the NREFT expectation. 

As shown by the dotted and dot-dashed lines in Fig. [Til we also observe a suppression of 
the meson loops in r( c — > XcO 71 " an d XcO Vc 7 ? which, however, is weaker as in the NREFT 
calculation. As already mentioned, this may be understood as part of the short-range physics 
is already mimicked by the form factors in the ELA. Such an uncertainty seems inevitable in 
the ELA and a measurement of the ratios of branching fractions would be less sensitive to 
it. The normalized widths of rj' c — > Xc0^° and XcO Vc 7 ? have a different ordering compared 
with the NREFT results, and the r]' c — > XcO^r appears to be one order of magnitude larger. 
This is because the mass of the r]' c is much closer to the DD* threshold than rj c , and the 
form factor parameterization adopted in the ELA makes the effect of the proximity more 
prominent than the NREFT. 

The other decay channels of interest are between the radial excitation S- and P-wave 
charmonia, such as h' c — > ip'Tr and x' c o ~^ 7 l'c 7T °- The NREFT predicts that the width of 
XcO ~^ ^c 71 " i s tree- level dominant, but the charmed meson loops may give a significant 
contribution, while the transition h' c — > ip'ir is dominated by the meson loops. We find that 
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the ELA result for h' c — > ip'n has a much larger normalized decay width considering loops 
only, which indicates the dominance of the loops for this transitions. It agrees qualitatively 
with the NREFT one quite well though has a larger value. The ELA result of the transition 
XcO ~^ Vc^ 1S smaller than that of h' c — > ip'ir while much larger than that of the other 
SP transitions. We consider this as a qualitative agreement with the NREFT, since in the 
NREFT the pattern is similar. 

This is consistent with our prospect that these two approaches would agree with each other 
given the dominance of the meson loops in the transitions. It should also be pointed out that 
in the ELA, singularities would appear when the mass of the initial charmonium is larger 
than the artificial threshold Ai + A2 introduced by the form factor. When the "threshold" is 
approaching, one typically gets larger results. The comparison of these two approaches turns 
out to be a useful tool for our understanding the properties of the ELA. 

5.3 Transitions between two P-wave charmonia 

The transitions between P-wave charmonia, especially for the Xcj decays, allow for an es- 
pecially sensitive test of the meson loop transitions, as the power counting analysis suggests 
that the meson loop contributions relative to the tree-level transitions will be enhanced sig- 
nificantly (c.f. tabled]). Thus, the comparison between NREFT and ELA will help restrict 
the possible parameter ranges in the ELA as well as study further the implications of the 
adopted form factor parameterization. 
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Figure 12: Parameter- free ratios among the decay widths of x' c i ~* Xcj^ considering only 
the meson loops in the ELA. 



In the case of x' c o ~> Xd^° and x'ci ~^ Xc0^°, our calculation shows that the normalized 
partial widths, T(x' c0 — > Xcl7T )/(<7i<7i)> and the mass-dependence of r(x' cl — > Xc0^°) / {9i9i ) 
have similar structures as those given by the NREFT, but the relative magnitude of the 
decay widths is different. The ELA predicts larger results for the widths of the transitions 
Xco ~~ ^ XciK and x'ci ~~ ^ Xcj^° than the corresponding values in the NREFT, except for the 
x'd ~~ ^ XcO 71 " which is opposite. 
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Figure 13: Ratios of the decay widths x' C 2 ~^ Xci^° and x' C 2 ~^ Xc2^° to F(x' c i ~~ ^ Xc2^°) in 
the ELA. 



A comparison is also made for these two methods in the predictions of the ratios, B^°q^ 
and R l °°i2, as defined by Eq. ([65]) . In Fig. [121 these two ratios are plotted in terms of a 
values. It is sensible to observe the stability of the ratios within the varying x' c \ mass region. 
In comparison with the NREFT results in Fig. ratio It^u* is consistent with the dashed 
line there even in magnitude, while ratio i?ioi2 exhibits an inverse relative magnitude between 
r(Xcl ~~ > XcO 71 " ) and T(x' c i ~~ Xc2^ )- Although there exist significant discrepancies between 
these two methods here, we emphasize that the flat behavior of the ratios again suggests some 
systematic model-independent features of the meson loop contributions. Especially, there are 
two kinks in i2]*o°2 i R both methods, which indicate the opening of the D°D*° and D + D*~ 
thresholds. In this sense, it is also listed as in qualitative agreement in Table [5j 

We further compare the ELA results for the P-wave decays, x' C 2 ~~ ^ Xci^° and x' C 2 ~^ Xc2^°, 
with the NREFT ones. For these two transitions, the contributing loops have Mi = Mi = D* . 
The mass of the initial state M Xc2 is smaller than the threshold 2Mo* , and much smaller than 
the artificial threshold 2(Md* + oAqcd)) and hence the unphysical singularities would have 
little effect. So one may expect the ELA would give similar results compared with the NREFT. 
The meson loops in the ELA give, 



r (x' c2 ' 
r (x' c2 



Xcivr°)ioo P 

* Xc27T°)loop 



OMgfg'" keV . 
: 0.19^f keV 



(75) 



with a = 1.64 as determined previously. They are consistent with those given by the NREFT 
in Eq. (I67p . In Fig. [T3l we present the ratios of the decay widths, -Rjjm anc ^ ^2212> as defined 
by Eq. (|68p . to compare with the solid and dashed lines in Fig. [8[ It is interesting to see the 
consistency of these two methods here in both mass evolution and magnitude, as expected. 

5.4 Summary of the ELA results 

As can be seen from Table [5l the NREFT and ELA results for all the transitions are in 
qualitative agreement. For both the SS and PP transitions the results for the two approaches 
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are in quantitative agreement as long as the decaying charmonium state is not sufficiently 
heavy to allow the form factor singularities to matter numerically. The observed level of 
agreement between these two approaches provides a further support for the power counting 
of the NREFT. It demonstrates nicely two important aspects of this effective field theory, 
namely, the applicability of the non-relativistic treatment in these charmonia transitions (note 
that the ELA formulations are relativistic) , and the dominance of the long-range pieces of the 
meson loops which is driven by the unitary cuts for the SS and PP transitions. At the same 
time, it provides additional confidence for applying the ELA to reactions where the NREFT 
is no longer applicable, such as in case of charmonium decays into light final states. 

The situation is different for the SP transitions. Most of these transitions are expected 
to be dominated by "short-range" physics, i.e. the tree-level contributions, in the NREFT. 
From the effective field theory point of view, the tree-level amplitudes serve two important 
duties: On the one hand, they provide a way to include the short-ranged quark dynamics. 
On the other hand, they may absorb the ultra-violet behavior of the meson loops and, if 
necessary, absorb their divergences. In the ELA, part of the short-range physics, which is 
completely omitted so far in the NREFT, is parameterized by the form factors. Thus, for the 
SP transitions we expect, and indeed find, that there are quantitative differences between 
these two approaches. However, it is important to note that the qualitative role of the loops 
is the same in both treatments. 

6 Summary 

In this paper, the effects of charmed meson loops in the transitions between two charmonia 
with the emission of one pion or eta have been systematically investigated. The power count- 
ing for the loop amplitudes within the framework of a non-relativistic effective field theory 
is given. The difference among the power counting estimates for the ratios of tree-level and 
loop contributions for the various transitions considered comes from the quantum number 
difference of the involved charmonia. It is found that the loops are enhanced in the transi- 
tions between two ,5- wave charmonia by a factor of 1/v and in the transitions between two 
P-wave charmonia by a factor of 1/v 3 . As a result, even if the A(3872) is a cc charmonium 
with J PC = 1 ++ , the dominant contribution to its decays into XcJ 71 " is given by the interme- 
diate charmed meson loops rather than the tree-level ones. For the transitions between one 
5-wave and one P-wave charmonia, because of the competition between 1/v 3 and q 2 /Mp one 
should analyze case by case. The loops are highly suppressed for the decays ?// — > h c ir° and 
7]' c — > XcO^ which have a small phase space, while the suppression for the decay x' c o ~~ ^ Vc^ 
is more moderate, and they are even enhanced in the decay h' c — > ip'ir . 

Among the loop dominated transitions, predictions for the ratios among the decay widths 
are given. In these ratios, the dependence of the unknown coupling constants is canceled. 
A detailed calculation in the framework of an effective Lagrangian approach is also given 
in comparison with the NREFT results. We find that the results from these two methods 
are qualitatively consistent with each other. Significant deviations appeared only when the 
results became sensitive to the particular form of the form factors used in the ELA, pointing 
at possible improvements for this approach, or when the considered transition is dominated 
by short-range physics. Although the lack of constraints on the structure of the counterterms 
would lead to uncertainties with the absolute magnitudes of the partial decay widths, the 
ratios between related channels are less sensitive to such uncertainties, and model-independent 
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aspects of the meson loop contributions can be highlighted. Those loop-dominant channels 
should be testable experimentally in the future, for instance, with Belle, PANDA, BES-III, 
LHC-b and Super B factory. 
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A Loop functions in dimensional regular izat ion 

Define the basic three-point scalar loop function in ci-dimension as 

1 



f drl 

I(q) = 1 J (2tt)<* (f-mj+ie) [( 



P-lf 



m; 



+ ie] [(I — q) 2 — m| + ie] 



(A.l) 



Non-relativistically, in the rest frame of the initial particle, i.e. = {M, 0}, it can be worked 
out analytically by performing the contour integration over I , and then integrating over the 
spatial components of the loop momentum using dimensional regularization, 
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(A.2) 



where fiij = rriimj/ (nii + rrij) are the reduced masses, 612 = mi + rri2 — M, 623 = rri2 + ms + 
q° - M, and 



f fJ-23\ ^2 r, , / r, , . M23 ^2 

a = ) q , c = 2^12612, c = 2^23023 H <? • 

V m-3 / 7773 



(A.3) 
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Since there is no pole for d = 4, in the last step, we have taken d = 4. 
We also need the vector and tensor loop integrals which are defined as 



2 - m\ + ie) [(P - I) 2 -m? 2 + ie] [(I - q) 2 - m§ + ie] ' 



(A.4) 



and 
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By using the technique of tensor reduction, we get 
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where the function B(c) is defined as 
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In the last step, the dimension d is taken to be 4. 



B Amplitudes for the charmonium transitions 

Here we give the amplitudes for the transitions between charmonia with the emission of one 
pion or rj discussed in the paper. Notice that the expressions for the amplitudes are only given 
for the charged charmed meson loops. The expressions for the neutral or strange charmed 
meson loops can be obtained easily by replacing the charged charmed meson masses by the 
the corresponding neutral or charmed ones. To distinguish loops with different mi, mi and 
m3, we write the loop functions as, for instance, I^\q, Ml, M2, M3) with Mi denoting the 
meson with the mass rrij in the following. 
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B.l Transitions between the S-wave charmonia 

• J/Y>7T (77) 

= »2|W2^' fc 9V(^)e fc ( J/^)^ 2 [2/} 2) (9, £> ± , J D ± , + 27p } (g, £*) 
+4/J 2) (g, D± , D** , + 24 2) (g, £>± , , - /« (g, D ± , D*± , D*± ) 

+4/J 2) (g, D**, D*±, £>*) + 2/<; 2) (g, D** J D ± ) - /«(g, D*±, D±) 
-2ij 2) (g, D* ± ,D ± ,D* ± ) - 21 {2) {q, D* ± ,D ± ,D* ± ) + /«(g, D**, D ± , 
- 10l[ 2) (g, D*± , D*± , D*±) - 2I {2) (g, D*± , D*± , D**) + (q, D*± , D** , D^) . 

(B.l) 

B.2 Transitions between the S- and P-wave charmonia 

-iM (V' -> /i c vr°) ± 
= -*2 J;W 2 {? • ^')? • W (g, £> ± , D ± , D*±) - iW (g, D* ± ,D ± ,D* ± ) 
+/ (1) (g, D ± ,D* ± ,D* ± ) - I« (g, D*±, D*±, D*±) 
+g 2 e{^') • e(h c ) \l (1) (q, D* ± ,D ± ,D* ± ) + /«(g, L>*± £>±) 
-/ (1) (g, D ± ,D* ± ,D* ± ) - 7«(g, } . (B.2) 
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J/V) ± 

*f;3i<?2 {? • e(h c )q- e(J/V0 [21^ (q, D* ± ,D ± ,D ± ) - I(q, D* ± ,D ± ,D^ 

-2lV> (g, D* ± ,D ± ,D* ± ) + I(q, D* ± ,D ± ,D* ± ) 
+2I<M (g, £)* ± ) D*± ) £) ± ) _ j( g> £)*±, D*±, D±) 

-2/W (g, £)*± ) £)*± ) £)*±) + £>*±, _D*±, £)*±) 

+2/W (g, D ± , D**, - I(q, D ± ,D* ± ,D* ± ) 
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-2/« (g, D** D**) + /(g, _D*±)j J . 



(B.4) 
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• XcO ->• ?? c vr° 
- %M (Xco -> 

B.3 Transitions between the P-wave charmonia 

• XcO -> Xcl^T 

- iM (x'co -> Xcivr ) 

• Xci -> Xcovr 

• Xcl -> XclTT 

- *M (Xci -> Xci^°) ± = -*2| ; Wie ij ' fc ?V(x' cl )e fc (xci)/(?,i5* ± ,I> ± ,D* ± ). (B.8) 

• Xcl -> Xc27T° 

- iM (x'ci -> Xc2vr°) ± = z2^| 5l < 7 ' 1 e^( Xc2 ) (7 %^( X ' cl )/( (? , J D ± , J D* ± , J D* ± ), (B.9) 

where £ u (xc2) is the polarization tensor of the Xc2, and it is traceless, i.e. £ n (xc2) = 0, 
and symmetric. 

• X c2 -> Xcivr 

- *M (x c2 -> X c ivr°) ± = -i2^| ffl < 7 / ie ^( Xc2 )^^( X ^ 1 )/(g, J D* ± ,Z)* ± , J D ± ). (B.10) 

• X c2 -> Xc2vr° 

-iM (X C 2 X c2 vr°) ± = -i4| 5l ^ e ^gV 7 ( Xc2 )^(Xc2)/(^ J D* ± , J D* ± , J D* ± ). (B.ll) 

• /l' c ->■ /l c 7T° 

-iM(^->/ lcir 0) ± = -il.g ig ' ie ^ q ^(h' c )e\h c ) 

x [J(g, D* ± ,D ± ,D* ± ) + J(g, D ± ,D* ± ,D* ± ) 

-I(q, D**, D**, D ± ) + J(g, D**, . (B.12) 
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C The Effective Lagrangian Approach (ELA) 

C.l Effective Lagrangians 

Based on the heavy quark symmetry \<59\ I7U| . the relevant effective Lagrangians used here 
are: 

£l = i^-TriP^H^Hu] +R.c, (C.13) 
£ 2 = i^Tr[R^H 2 ^d„H^ +H.c, (C.14) 

where the spin multiplets for these four P-wave and two S*-wave charmonium states are 
expressed as 

2 

(C.15) 



±iy V *^-V il5 ], (C.17) 
ff* = 1©,*% - A75] f ^) , (C18) 



Re = [-^)(^-Vcj 5 ){-^)- (C.16) 

The charmed and anti-charmed meson triplets read 

li 



where T> and T>* denote the pseudoscalar and vector charmed meson fields, respectively, i.e. 
£>(*) = (_ D o(*) j D +(*) } D+W). These La grangians can be reduced to the two-component ones 
used in the NREFT in Section El 

Consequently, the Lagrangian for S- wave J/tjj and r/ c is 

+9 Vc D* D D^(d^ c D - rjcd^D) + ig VcD *D*e^d^DlD^d^ (C.19) 

and the Lagrangian for P-wave h c and Xcj is 

£p = i9x*D*D*x"lKD$ + g Xcl D*DXd^D + ig XcoD * D *D*»D*J + ig XcoDD DD 

+g hc D* D D^h^ + ig hcD * D *e^ aP Dld a h cu Dl + H.c, (C.20) 

where the couplings are 



g-^DD = g2^niipm D , g^D*D* = -g2^m 1 i > m D *, g^D*D = g2\l ^ ■ 



g Vc D*D = g2^rn Vc m D , g Vc D*D* = 52 
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gh c D*D = -giy/m hc m D m D *, gh c D*D* = 9i\ 



m2 p* 
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(C.21) 



In the following, we present the transition amplitudes for the intermediate meson loops 
listed in Table[3]in the framework of the ELA. Some conventions should be clarified in advance. 
Notice that the expressions are similar for the charged, neutral and charmed strange mesons 
except that different charmed meson masses are applied. We thus only present the amplitudes 
for those charged charmed meson loops. G1G2G3 is the product of vertex couplings for each 
loop, and the explicit expression can be found in Eq. (ICl2ll . T = Ui^ii^Qi) is the 
form factor, [gi, 03, 02] are the four-vector momenta for the intermediate mesons [Ml, M3, 
M2], respectively, pi (pf), s% (£/), and fa (cj)f) are the initial (final) charmonium four-vector 
momentum, polarization vector, and polarization tensor, respectively. P2 is the four-vector 
momentum of 7T° or rj. 

C.2 Transitions between the S'-wave charmonia 
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C.3 Transitions between the S- and P-wave charmonia 
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C.4 Transitions between the F-wave charmonia 

• XcO -> XclTT 
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